Self-localization of magnon Bose-Einstein condensates on the ground and excited 

levels: from harmonic trap to a box 
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Long-lived coherent spin precession of 3 He-B at low temperatures around 0.2 T c is a manifestation 
of Bose-Einstein condensation of spin-wave excitations or magnons in a magnetic trap which is 
formed by the order-parameter texture and can be manipulated experimentally. When the number 
of magnons increases, the orbital texture reorients under the influence of the spin-orbit interaction 
and the profile of the trap gradually changes from harmonic to a square well, with walls almost 
impenetrable to magnons. This is the first experimental example of Bose condensation in a box. By 
selective rf pumping the trap can be populated with a ground-state condensate or one at any of the 
excited energy levels. In the latter case the ground state is simultaneously populated by relaxation 
from the exited level, forming a system of two coexisting condensates. 
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During the last few years increasing efforts have been 
invested in the investigation of Bose-Einstein condensa- 
tion (BEC) of non-conserved bosons, such as magnons 
[H, [![, photons Q, and exciton-polaritons Q. BEC 
of quasiparticles and other particle-like excitations is a 
special case since in thermal equilibrium their chemi- 
cal potential vanishes. Formally BEC requires conser- 
vation of the particle number, but condensation can still 
be extended to systems with weakly violated conserva- 
tion if a dynamic steady state is created. The loss of 
(quasi)particles owing to their decay can be compensated 
by pumping and thus, for sufficiently long-lived excita- 
tions, the non-zero chemical potential is well defined and 
condensation becomes possible. Magnons in supcrfluid 
3 Hc-B satisfy this condition and condensation there is 
observed as spontaneous long-lived coherent precession 
of spins which is accompanied by various phenomena of 
spin superfluidity [f|. Different types of magnon BEC 
have been identified in 3 He-B, starting from the so-called 
"Homogeneously Precessing Domain" (HPD), which rep- 
resents a bulk condensate state. 

In 3 He-B at low temperatures the HPD state becomes 
unstable owing to parametric creation of spin waves, but 
another type of coherent long-lived NMR signal with sev- 
eral orders of magnitude smaller amplitude was discov- 
ered in this regime @ . The mode was ascribed to magnon 
condensation in a magnetic trap formed in a weakly in- 
homogeneous order parameter distribution or texture [l[ , 
but the piecemeal information obtained with pulsed and 
cw NMR measurement has been confusing. Here we re- 
port the first measurements with full experimental con- 
trol of the order parameter texture which produces the 
3D trap. We find that when the number of magnons in- 
creases the profile of the trap changes from harmonic to 
a box. The pressure of the multi-magnon wave function 



opens a "cavity" in a way similar to the electron bubble 
in liquid helium. In quantum field theory such self local- 
ization of a bosonic field is known as a Q-ball 0. This 
texture-free "cavity" can be filled by a magnon conden- 
sate on any energy level of the trap. For cold atoms in 
an optical trap, the formation of a non-ground-state con- 
densate has been discussed Q, but not yet realized. Of 
great practical importance is the fact that these magnon 
condensates can be used to probe the quantum vacuum 
state of 3 He-B in the limit T — > [§], where most con- 
ventional measuring signals become insensitive. 

Magnon condensation:-There are two approaches 
to the thermodynamics of atomic systems: one can fix 
the particle number J\f or the chemical potential fx. For 
magnon condensation, this corresponds to different ex- 
perimental situations: to pulsed or continuous wave (cw) 
NMR, respectively. In free precession after the tipping 
pulse, the number of magnons pumped into the trap is 
conserved (if losses arc neglected). This corresponds to 
fixed Af, when the system itself chooses the global fre- 
quency of coherent precession (= the magnon chemical 
potential /x Q). The opposite case is cw NMR, when a 
small rf field is continuously applied to compensate for 
the losses. The frequency of precession uj is then that of 
the rf field w r f , the chemical potential is fj, = uj = cj r f , and 
the number of magnons adjusts itself to this frequency, 
to match the resonance condition. 

A cylindrically symmetric trap for magnons is schemat- 
ically shown in Fig. [TJ It is realized in a long cylindrical 
sample container with radius R s ~ 3 mm in an axially 
oriented magnetic field (the experimental setup is de- 
scribed in Ref. [l(|). The axial confinement potential 
Uh(z) = ujl(z), where wl(z) = jH(z) is the local Lar- 
mor frequency, is produced by a small pinch coil, which 
creates a shallow minimum in the magnetic field. The 
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FIG. 1: A sketch of the trapping potential Eq. @ which is 
formed in the cylindrically symmetric "flare-out" texture of 
the orbital anisotropy axis L (thin lines) in a shallow min- 
imum of the vertical magnetic field H (right). The arrows 
represent the magnetization M, which precesses coherently 
with constant phase angle in the condensate droplet (dark 
blue), in spite of the inhomogeneity in the texture and in the 
magnetic field. 



radial confinement comes from the spin-orbit interaction 
with the texture of the orbital angular momentum L of 
Cooper pairs 
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where VP is the wave function of the magnon condensate, 
51l is the Lcggctt frequency characterizing the strength 
of the spin-orbit coupling, and /3l the polar angle of L. 
The trapping potential is then given as 



U(T) = U l \(z) + U ± {r)=w L (z) + 



5w L 2 
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On the side wall of the cylinder L is oriented normal to 
the wall, while on the cylinder axis L || H. Close to the 
axis /?l remains small and varies linearly with distance 
r. Here the potential ?7(r) reduces to that of a harmonic 
trap , as used for the confinement of dilute Bose gases 
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U(r) = U(0) 



m M 
2 



2 2 
U1ZZ 



2 2\ 



(3) 



where mM is the magnon mass. The low-amplitude 
standing spin waves have the conventional spectrum 



w L (0) + u r (m + 1) + w 2 (n + 1/2) . 



(4) 



where wl(0) = ul(z = 0) is the Larmor frequency at the 
bottom of the well, which corresponds to the center of 
the trap in Fig. [TJ The axial oscillator frequency u> z is 
adjusted by changing the current in the pinch coil, while 
the radial frequency u r can be controlled by rotating the 
sample, since adding vortex-free supcrfluid flow or rec- 
tilinear vortex lines modifies the flare-out texture. Let 
us consider the condensates which form when we start 
filling magnons to one of the levels (m,n) in Eq. (0]). 

Ground-state condensate:-When the number of 
magnons N in the ground state (0,0) increases, they ex- 
ert an orienting effect on the L texture via the spin-orbit 
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FIG. 2: Calculations of multi-magnon bubbles with the con- 
densate in the ground state, n — m — 0. The deflection angles 
of the magnetization, /3m, and of the textural anisotropy axis, 
/3l , are plotted as a function of radius for different condensate 
populations. The population increases from bottom to top in 
the upper plot and from left to right in the lower plot. The red 
curves correspond to the experimental regime in Fig. [3] (start- 
ing with Mx/Mhpd = 5-1CP 4 and then continuing from 10 -3 
to 5 • 1CP 3 with a step of 10 -3 ). The blue curves illustrate ex- 
trapolations to the asymptotic regime / — > /l(0) (Mx/A/hpd 
varies from 1CT 2 to 5 • 1CT 2 with a step of 5 ■ 10~ 3 ). When 
the magnon occupation increases, the magnon wave function 
suppresses the orbital texture /3l and the potential well trans- 
forms towards a box with impenetrable walls. Fits to the 
wave function of the condensate in a box are shown with bro- 
ken lines in the upper plot. The effective radius of the box 
obtained from such fits is shown in the inset as a function of 
the magnon occupation number A/". The slope from Eq. 
with p — 1/5 is shown for comparison. 



interaction in Eq. ([T]), which favors L || H. As a result 
at large Af the harmonic trap transforms to a box with 
/3l ~ within which magnons are localized. This effect is 
demonstrated in Fig. [2] with self-consistent calculations of 
the texture and of the magnon condensate wave function 
"J in the axially symmetric and z-homogeneous geome- 
try. Rapid flattening of the texture and self localization 
of the wave function is seen to result. At large M the 
radius of localization approaches the asymptote 



(5) 



where a r is the harmonic oscillator length in the original 
radial trap (at J\f <C A/" c ), M c is the characteristic num- 
ber at which the scaling starts, and p « 0.2. As in the 
case of the electron bubble, R is determined by a balance 
between the magnon zero-point energy and the surface 
energy of the condensate bubble. For the 2D radial tex- 
ture the total energy is: 

The first term on the rhs is the kinetic energy of Af 
magnons in a cylindrical box, where A m is the m-th root 
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FIG. 3: Formation of the magnon condensate in the ground 
state of the trap in Fig.[T]in cw NMR measurement. The con- 
densate magnetization M± precessing in the transverse plane 
is plotted on the vertical axis, normalized to that when the ho- 
mogeneously precessing domain (HPD) fills the volume within 
the detector coil [l2l ] . The arrows indicate the sweep direction 
of the applied rf frequency / = w/(2tt). M± grows when the 
frequency is swept down. Only a tiny response is obtained 
on sweeping in the opposite direction. During the downward 
frequency sweep the condensate is destroyed (vertical lines), 
when energy dissipation exceeds the rf pumping. This point 
depends on the applied rf excitation amplitude, marked at 
each vertical line. The lower green line represents the result 
of calculations from Fig. [2] These calculations have no fitting 
parameters: the vertical difference from the measurements 
can be attributed to the experimental uncertainty in deter- 
mining the normalization for M±. (Insert) The experimental 
curve measured with the largest excitation and compared to 
the numerical curve from the main panel, replotted with log- 
arithmic axes to demonstrate the asymptotic limit for large 
magnon numbers: M± oc (/ — /l(0)) -1 , which corresponds 
to the condensate in a box. 



of the Bessel function. The second term is the surface 
energy with the surface tension a which depends on R 
due to the flexibility of the texture. Our numerical sim- 
ulations give cr(R) oc i? 2 , and minimization of Eq. Q 
with respect to R gives for the box radius Eq. ([5]) with 
pral/5. 

Incidentally, for an atomic condensate with repulsive 
inter-particle interactions the exponent in Eq. ([5]) is also 



p = 1/5 in the Thomas- Fermi limit [1 1| . However, owing 
to different scenarios of condensate formation the depen- 
dence of the frequency shift on TV" differs from the behav- 
ior of the analogous quantity in an atomic condensate, 
i.e., the chemical potential fi(Af) is of the form 

w-Wl(O) ~ w r (M/M c )~ 2/5 , magnon BEC, (7) 
fi-U(0) « w r (M/M c f' 5 , atomic BEC. (8) 

In contrast to condensates of ultra-cold atoms, in the 
magnon condensate dw/dN < 0. When the magnon con- 
densate is growing, its frequency w decreases, approach- 
ing the Larmor frequency wl asymptotically. This dctcr- 
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FIG. 4: Formation of magnon condensates at different ex- 
cited levels (m,n) of the magnetic trap in Fig. [TJ When the 
frequency of the rf excitation field is swept down and crosses 
one of the levels (m,n) at / m n, the magnon condensate starts 
to grow according to Eq. Q. Only modes with even m and 
n are excited. The three examples of measured spectra are 
for different values of rotation velocity Q in the vortex-free 
state and for different depths AH/H of the axial field mini- 
mum in the trap. In the topmost case the depth of the field 
minimum is twice that of the two bottom examples. As a 
result, the smaller axial level spacing increases, as the axial 
oscillator frequency oj z /2n : 22 — > 27 Hz, but the larger ra- 
dial level spacing remains unchanged. Similarly, when 51 is 
increased from 0.6rad/s (bottom spectrum) to 0.8 rad/s (two 
top most spectra), the radial level separation increases, as 
uj r /2TT : 170 — > 240 Hz. The equidistant vertical tick marks 
refer to the level positions in the absence of magnons and have 
been fitted to the harmonic trap relation in Eq. Q. 



mines the way in which the magnon condensate is grown 
in a cw NMR experiment (Fig. [3]). Magnons are cre- 
ated when the frequency uj of the applied rf field is swept 
down and crosses the ground state level woo- With in- 
creasing N ', the potential well becomes wider radially 
and the energy of the trapped state decreases, finally ap- 
proaching the scaling regime. Assuming that a(R) cx R 2 
and taking into account that M = J d 2 r\^>\ 2 and thus 
|\Jr| ~ A/" 1 / 2 / R, one obtains for the transverse magnetiza- 
tion Mi cx fd 2 r\V\ (xN 1 ' 2 R oc R 1 / 2 oc (/ - /l(0))" 7/4 . 
This agrees with the exponent —1.75 obtained in numer- 
ical simulations and in the experiment (Fig.[3j insert). 

Non-ground-state condensates:-In similar manner 
the excited states (m, n) can be populated and the entire 
spectrum scanned: the condensate in the state (m, n) 
starts to grow when the rf frequency is swept down and 
crosses the level w mn (0) from above, since with increasing 
Af the frequency decreases. In the limit of large magnon 



4 



10 20 Time, s 




0.2 0.4 0.6 0.8 1 1.2 



/-/ L ,kHz 

FIG. 5: Decay of the magnon condensate, after rf pumping is 
switched off at t = 0. The amplitude of the Fourier transform 
of the signal from the NMR pick-up coil is shown in the main 
panel. Magnons are pumped to the (2,0) level at t < 0, but the 
ground state is simultaneously populated owing to the decay 
of magnons from the excited state. At t > both states decay 
and the frequency of precession increases as the trap responds 
to the decreasing magnon population (Fig. |3J. The insert 
shows the amplitudes of the two peaks and fits to exponential 
decay with time constant r. The measurements have been 
performed at 0.14T C in equilibrium rotation at 0.9rad/s with 
rectilinear vortices. 

occupation M in an excited state, the resulting Q-ball 
is analogous to a bubble in liquid He with an excited 
electron. 

Three examples of level spectra are shown in Fig. [4] 
The figure demonstrates how the textural trap is con- 
trolled by rotation in the vortex-free Landau state, where 
the velocity of the superfluid fraction is v s = 0. Here 
the azimuthally flowing superfluid counterflow velocity 
v n — v s = |S~2 x r| modifies the L-texture, making the trap 
steeper in the radial direction. This increases the oscil- 
lator frequency uj r and the spectral distance between the 
radial modes (m). Similarly, on increasing the depth of 
the minimum in the field H the distance between the 
axial modes (n) increases. 

The condensates can be maintained in steady state by 
continuous pumping at ui mn (J\f). After switching off the 
pumping, the condensate is manifested as long-lived ring- 
ing of the free induction signal. In Fig. [3] the decay of 
the induction signal has been recorded, showing two co- 
existing condensates: at the excited (2,0) level, where the 
magnons were initially pumped, and at the ground level 
(0,0), which is filled by relaxation from the excited state. 
This quantum relaxation process is similar to the forma- 
tion of the magnon condensate with incoherent pumping 
and explains the off-resonance excitation of the ground 
state population observed in Ref . [l3| • During the decay 
the magnon population follows closely the trajectory for 
the reverse process of Fig. [3] With decreasing temper- 



ature the relaxation rate in the ground state decreases 
rapidly life times ~ 15min have been rep orted from 
observations at the lowest temperatures [14| . Our mea- 
surements in Fig. [5] on the decay from the excited state 
(2,0) show faster relaxation. However, the measured re- 
laxation time is much longer than the dephasing time 
of the linear NMR response (about 10 ms in the same 
conditions). This coherence of the precession is the most 
striking experimental signature of the Bose-Einstein con- 
densate, not only in the ground state, but also in the 
excited states of the trap. 

Conclusions:-Wc have demonstrated the forma- 
tion of coherently precessing magnon condensates in a 
magneto-textural trap with experimentally controllable 
potential. As distinct from the traps of ultra-cold atoms, 
here the trap transforms with increasing magnon num- 
ber from a harmonic well to a cylindrical box. Different 
excitation levels can be selectively populated with con- 
densates. These provide new opportunities for measure- 
ments in the T — > limit Q, particularly in terms of 
their relaxation properties. An urgent task is to deter- 
mine whether the relaxation rate will reveal new infor- 
mation about surface fljj and vortex-core bound states 
[lit which are Majorana-fermion-like zero-energy modes 
of the topological insulator superfluid 3 He-B. So far large 
temperature-independent surface relaxation has been re- 
ported when the magnetic trap borders to a boundary 
[1 7I | , while we have preliminary observations of enhanced 
relaxation with increasing vortex number within the con- 
densate. The proper understanding of these effects re- 
mains a task for the future. 
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